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ABSTRACT 

We present the first sample of 31-GHz selected sources to flux levels of 1 mJy. From late 2005 to 
mid 2007, the Sunyaev-Zel'dovich Array (SZA) observed 7.7 square degrees of the sky at 31 GHz to a 
median rms of 0.18 mJy/beam. We identify 209 sources at greater than 5a significance in the 31 GHz 
maps, ranging in flux from 0.7 mJy to ~ 200 mJy. Archival NVSS data at 1.4 GHz and observations at 
5 GHz with the Very Large Array are used to characterize the sources. We determine the maximum- 
likelihood integrated source count to be N{>S) = {27.2 ± ^.5)deg~^ ^{Smjy)~^ ''^^"''^ over the 
flux range 0.7 — 15 mJy. This result is significantly higher than predictions based on 1.4-GHz selected 
samples, a discrepancy which can be explained by a small shift in the spectral index distribution for 
faint 1.4 GHz sources. Prom comparison with previous measurements of sources within the central 
arcminute of massive clusters, we derive an overdensity of 6.8 ± 4.4, relative to field sources. 
Subject headings: techniques: interferometric, catalogs, surveys, cosmology: cosmic microwave back- 
ground, cosmology: observations, radio continuum: general 



1. INTRODUCTION 

Until recently, prohibitive integration times have lim- 
ited deep surveys at high radio frequency to small areas, 
or high flux cutoffs. As a result, estimates of the proper- 
ties of compact extragalactic radio sources at frequencies 
above ~ 15 GHz have come primarily from observations 
at low- frequencies (in the < 10 GHz range), where the 
bulk of the sources are brighter, extrapolated to higher 
frequencies by targeted followup campaigns (de Zotti 
et al. 2005). While an understanding of high-frequency 
source populations is interesting in its own right, it is 
also important for the current generation of CMB and 
SZ experiments, many of which operate at frequencies 
> 30 GHz, where the source population is poorly char- 
acterized. 

In recent years the development of broad-band corre- 
lators has made deep surveys of significant areas of sky 
possible with interferometric arrays operating at 31 GHz. 
In this paper, we report results of a 7.7 square degree, 
31-GHz sky survey with the Sunyaev-Zel'dovich Array 
(SZA). Although the primary goal of the SZA was to 
measure CMB anisotropy and to search for galaxy clus- 
ters via their Sunyaev-Zel'dovich (SZ) effect, it has also 
yielded the first catalog of ^mJy sources selected at 
31 GHz. This paper focuses on the results of the SZA 
survey observations as they pertain to the population of 
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high-frequency (^ 31 GHz) selected compact sources. 

Several experiments, such as WMAP (Bennett et al. 
2003), DASI (Kovac et al. 2002), CBI (Mason et al. 
2003), and the VSA (Cleary et al. 2005), have character- 
ized sources at 31 GHz brighter than 5 mJy. The present 
work extends this characterization to the fainter 31 GHz 
sources. As these faint sources are at or near the noise 
level of high-frequency CMB experiments, they repre- 
sent a serious contaminant which must be carefully mod- 
eled and statistically accounted for. The OVRO/BIMA 
SZ group (Coble et al. 2007) have surveyed ~ 29-GHz 
sources associated with massive clusters at the mJy 
level; the SZA survey allows a comparison of this highly- 
selected population with the bulk properties of field 
sources at 31 GHz. More recently, the CBI collaboration 
has followed up over 3000 1.4 GHz-selected sources at 
31 GHz (Mason et al. 2009) and used these observations 
to predict the source population at 31 GHz. Comparison 
of this prediction with the SZA measurement at 31 GHz 
provides a direct test of the assumptions that underlie 
this extrapolation. 

The paper is organized as follows: in §2 we present 
a description of the instrument and of the SZA observa- 
tions. In §3 we describe follow-up observations performed 
with the Very Large Array (VLA), while §4 details the 
algorithm used to extract source fluxes from the SZA 
survey data. The characteristics of the 31 GHz-selected 
sample of sources are presented in §5, followed by a dis- 
cussion in §6. Conclusions are presented in §7. 

2. SUNYAEV-ZEL'DOVICH ARRAY OBSERVATIONS 

2.1. The Sunyaev-Zel'dovich Array 

The Sunyaev-Zel'dovich Array is an interferometer de- 
signed specifically for detecting and imaging the SZ effect 
in galaxy clusters, and is located at the Owens Valley Ra- 
dio Observatory (OVRO). The SZA is equipped with an 
8-GHz wideband correlator and sensitive 26-36 GHz and 
85-115 GHz receivers. In this paper, we present results 
only from 26-36 GHz (hereafter 31 GHz) SZA observa- 



2 



tions. 

The SZA consists of eight 3.5-meter antennas. For the 
observations presented here, six were arranged in a close- 
packed configuration (yielding high brightness sensitivity 
on angular scales typical of clusters of galaxies) , and two 
outlier antennas provided long baselines for sensitivity to 
compact objects. The SZA can therefore be thought of 
as two complementary interferometers: one with a typ- 
ical resolution of a few arcminutes (short antenna sepa- 
rations, or baselines), and one with resolution of about 
23 arcseconds (long baselines). For a more detailed dis- 
cussion of array layout and corresponding resolution, see 
Muchovej et al. (2007). 

In the limit where sky curvature is negligible over the 
instrument's field of view, the response of an interferom- 
eter on a single baseline, known as a visibility, can be 
approximated by: 

V{u, v)^ A{1, m)I{l, m)e-2"W+'""}d? dm, (1) 

where u and v are the baseline lengths projected onto 
the sky, I and m are direction cosines measured with re- 
spect to the (u, v) axes, A(l, m) is the normalized antenna 
beam pattern, and 1(1, m) is the sky intensity distribu- 
tion. Eq. 1 is a two-dimensional Fourier transform, the 
inverse of which is the image of the source intensity mul- 
tiplied by the primary beam pattern, known as a dirty 
map Id: 

loil, "ni) = A{1, m)I{l, to) 

p p+ca 

= V{u,v)e'^'''^''^+''"'Uudv. (2) 
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In practice, an interferometer measures discrete Fourier 
modes, and structure in the dirty map is convolved with 
a function which reflects this incomplete Fourier-space 
sampling. This function, called the synthesized beam, is 
equivalent to the point-spread function for the interfer- 
ometer. 

It is clear from Eq. 2 that the field of view for each 
dirty map is limited to the size of the primary beam, 
A{l,m), namely - 11.0' (FWHM) for the SZA at band 
center (30.938 GHz). To image larger areas, we use lin- 
ear mosaicking to stitch together nearby pointings which 
cover the region of interest (e.g., Pearson et al. 2003). 
Note that although large areas of sky can be surveyed in 
this manner, the resulting mosaics contain no informa- 
tion on size scales larger than that probed by the shortest 
baseline. 

The results presented in this paper are derived from 
two large SZA projects: one to measure primary CMB 
anisotropy (Sharp et al. 2009), and the other to survey 
for galaxy clusters via their SZ effect. Observations were 
conducted differently for these two projects, as detailed 
in §2.3. 

2.2. Field Selection 

The SZA survey fields were selected to lie far from 
the plane of the Galaxy and to transit at high elevation 
at the OVRO site, minimizing atmospheric noise while 
optimizing the imaging capabilities of the array. Fields 
were spaced equally in Right Ascension (RA) to permit 
continuous observation. These constraints led to the se- 
lection of four regions ranging in declination from 25° 



to 35°. Figure 1 depicts the approximate locations of 
these four fields overlaid on the IRAS 100 dust maps 
(Clegg 1980). 




Fig. 1. — IRAS 100 fim dust map with overlay of the SZA field 
locations. 



2.3. Observation Strategy 
2.3.1. Survey Mode 

Each of the four fields is split into 16 rows of 16 point- 
ings. The pointings are equally spaced by 6.6' along 
great circles in the RA direction, and each row is equally 
spaced by 2.9' in the DEC direction. Subsequent rows 
are offset from one another so that that the first pointing 
in each row is shifted by 3.3' in the RA direction relative 
to the previous row. This means that for a single field we 
observe an area that spans roughly 2 degrees in the RA 
direction and 1 degree in the DEC direction (see Figure 
2). 
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Fig. 2. — Mosaic pointing locations for a given SZA survey field. 
The fields are divided into 16 rows of 16 columns, with the pointings 
in each row separated by 6.6' and each row offset from each other 
by 2.9'. This leads to each field being roughly 2 degrees by 1 degree 
in area. In a single track the SZA observed four pointings within a 
given row. For example, pointings in the first and nintli, followed 
by pointings in the second and tenth columns. 



For each of the survey fields, data were taken daily in 
6 hour tracks. In a single track, we observed two stag- 
gered pairs of pointings, all within a single row. These 
observations were performed in a manner that permits 
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ground subtraction from consecutive pointings in a pair 
(although ground subtraction was not used in the analy- 
sis presented here). Each track results in roughly 1 hour 
of observation on each of the four pointings, with very 
nearly the same Fourier sampling for pairs of pointings. 
A second track is run at a later date, with the order of 
the pairs reversed, to ensure that the Fourier sampling 
for all four pointings is comparable. In Figure 2 we show 
the position of the pointings in each field, and indicate 
how the pointings were observed in a given track. 

For each set of four pointings, this sequence is repeated 
three times over the span of roughly one year, so that 
each pointing is observed in six total tracks, translating 
to roughly 6 hours of observation per pointing over the 
duration of the survey. 

2.3.2. CMB anisotropy Observations 

In addition to the survey observations, data were sep- 
arately taken to measure the anisotropy in the CMB. 
These consisted of observations of 44 distinct pointings, 
each separated by one degree, which were not mosaicked, 
but analyzed individually. Of these 44 pointings, 11 over- 
lap with pointings in the survey fields described above, 
and the rest are within a two degree radius of the center 
of the four survey fields. Where they overlap, the anal- 
ysis in this paper uses the survey data only. The track 
structure in the anisotropy observations is similar to that 
in the survey analysis; see Sharp et al. (2009) for further 
details. 

2.4. Observations 

Table 1 presents details of the mosaicked SZA survey 
observations (see Sharp et al. (2009) for the equivalent 
information on the CMB anisotropy observations). The 
second and third columns show the approximate pointing 
center of each 16-row field. We also present the bandpass 
and gain calibrators in the next two columns, with their 
fluxes as measured by the SZA. In the fifth column we 
give the time range over which observations were taken, 
with the caveat that observations were not performed ev- 
ery day during that time span. The penultimate column 
lists the total unfiagged integration time for data used 
in the analysis, and the final column gives the number 
of rows observed in each field. To ensure uniform cov- 
erage of all fields, tracks were repeated when necessary. 
Note that the full 16 rows were not observed for all fields, 
due to maintenance operations, instrumental character- 
ization, and RFI monitoring. For the first 8 months of 
observations, the SZA4 field was used for the dedicated 
CMB anisotropy measurements described above. As a 
result, only 7 rows in the SZA4 field were completed in 
survey mode. 

The data in the SZA survey correspond to 1493 tracks 
taken between November 13, 2005 and July 25, 2007. 

The data in the CMB anisotropy measurements corre- 
spond to an additional 414 tracks taken between Novem- 
ber 12, 2005 and October 25, 2007. The analysis in this 
paper refers to the full 1907 tracks taken in both observ- 
ing modes. 

Data for an individual track were calibrated using a 
suite of MATLAB^ routines, which constitute a complete 

9 The Mathworks, Version 7.0.4 (R14), 

http : //www . mathworks . com/products/matlab 



pipeline for flagging, calibrating, and reducing visibility 
data (Muchovej et al. 2007). Although the data were re- 
duced exactly as described in that paper, data collection 
differed in a few key ways: four distinct pointings were 
observed before observing a calibrator, and system tem- 
perature measurements were performed every eight min- 
utes. The absolute flux calibration is referenced to Mars, 
assuming the Rudy (1987) temperature model, and is es- 
timated to be accurate to better than 10%. Typical sys- 
tem temperatures measured throughout the survey were 
in the range 40-50 K. Flagging of the data as described 
in Muchovej et al. (2007) resulted in a loss of roughly 
15% of the data. At the end of a single 6-hour track, we 
achieved a noise level of approximately 1.5 mJy/beam in 
each pointing of the short and long baseline maps. 

2.5. Resulting Mosaics 

Once data on all pointings in a given fleld are reduced, 
we construct a linear mosaic of the field on a regular 
grid of 3.3" resolution. This sc;ale is much less than the 
requirement for Nyquist sampling of the data, ttft — ) 
where Dmax is the longest baseline, and ensures that the 
number of pixels over 2 degrees is a power of 2, con- 
venient for fast inversion of the Fourier data via FFT. 
The maps are a result of combining the data across our 
8 GHz of bandwidth, so that they approximate the sky 
at the central observing frequency. 30.938 GHz. The pri- 
mary beam is calculated from the Fourier transform of 
the aperture illumination of each telescope at the central 
observing frequency, modelled as a Gaussian with a cen- 
tral obscuration corresponding to the secondary mirror. 

Unlike the maps of individual pointings, the mosaicked 
maps are an estimate of the true sky signal at each point 
in the map; that is, the taper of the primary beam has 
been divided out. Due to the overlap of neighboring 
pointings, the effective noise is approximately uniform 
in the interior of the mosaics, but increases significantly 
towards the edge of the mosaicked images. We limit the 
survey area by applying an edge cutoff in our mosaicked 
maps where the effective noise is > 0.75 mJy/beam (cor- 
responding roughly to the one-third power point of the 
beam, given the noise in a single pointing). 

In Table 2 we show the noise properties of the observed 
fields. We present the minimum and median noise (in 
mJy/beam) for mosaic maps made with long baselines 
only, short baselines only, and with the combination of 
the two. The median noise is calculated only in the region 
within which the noise is less than the 0.75 mJy/beam 
cutoff. The last column indicates the total area covered 
in each field. That the minimum and median pixel noise 
values are similar is an indication of the uniformity of 
the coverage in the survey fields. This is not the case in 
the CMB anisotropy fields as they consisted of discrete 
pointings which do not overlap on the sky. 

3. VLA OBSERVATIONS 

As described above, the primary goals of the SZA 
project were a small-scale CMB anisotropy measurement 
and an SZ survey for clusters of galaxies, both of which 
require an accurate accounting of foreground emission. 
Although the long and short baseline data provide some 
intrinsic ability to discriminate compact objects, as dis- 
cussed in §2, high-sensitivity follow-up observations of 
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TABLE 1 
Survey Observations 



Field Name Field Center (J2000) Calibrators Dates Integration Rows 

a 5 Bandpass (Jy)^ Gain (Jy)* of Observations Time (hrs) Covered 



SZA4 02''15'"38'.3 32°08'21" J2253+161 (11.6) J0237+288 (2.9) 07/11/2006 to 07/25/2007 687 7 

DLS 09''19'"40=.0 30°0l'26" J0319+415 (11.0) J0854+201 (5.4) 11/18/2005 to 07/06/2007 1054 14 

NDWFS 14''30'"08'.0 35°08'34" J1229+020 (25.3) J1331+305 (2.1) 11/19/2005 to 07/23/2007 1000 14 

SZA3 2l'*30'"07^0 25°0l'26" J1642+398 ( 5.5) J2139+143 (1.4) 11/13/2005 to 07/25/2007 1245 16 



^Fluxes obtained from 31 GHz SZA observations of sources on April 16, 2006. 



the SZA fields were obtained with the VLA^° to faciH- 
tate source removal. 

VLA data at 1.4 GHz are publicly available from the 
NRAO VLA Sky Survey (NVSS; Condon ct al. (1998)) 
for all of the SZA fields, but are limited by the relatively 
coarse resolution of the NVSS, 45" FWHM, and high 
detection threshold of 2.5 mJy. The finer resolution (5") 
and deeper sensitivity {rms of 0.15 mJy) obtained with 
the FIRST survey (White ct al. 1997) arc better suited 
for this analysis, but data are only available on half of 
the fields (namely the DLS and NDWFS fields). 

To complement the NVSS and FIRST observations, we 
obtained high-sensitivity VLA observations at 5 GHz. 
Approximately 116 hours of observation were required 
to cover all four of the SZA fields, in four time blocks 
between 24 Feb 2007 and 15 April 2007. Data were taken 
at a center frequency of 4.86 GHz with the VLA D-array 
configuration, by mosaicking 180 pointings in each field. 
These pointings were arranged in 9 rows of 20 pointings 
in a hexagonal pattern, equally spaced by 6 arcminutes. 

The pointings were mosaicked, after CLEANing, us- 
ing the AIPS package FLATN, with the noise calculated 
using RMSD. Sources were extracted using the AIPS 
SAD algorithm, which itcrativcly removes the bright- 
est point in a mosaic using an elliptical Gaussian model. 
The rms achieved on the VLA 5 GHz observations was 
roughly 70 /zJy/beam, resulting in 859 sources from all 
four fields, down to a 5a depth of ^ 0.33 mJy. For the 
fields observed in CMB anisotropy mode, VLA follow-up 
at 8 GHz was obtained. These observations are described 
in Sharp et al. (2009). In this analysis, we use the VLA 
data to improve extraction of compact sources from our 
data as described in §4. 

4. SOURCE EXTRACTION FROM THE SZA SURVEY 

4.1. Overview 

Source identification begins in the image plane, with 
inspection of the combined (short and long baseline) sig- 
nificance (snr) maps for the brightest pixel with signifi- 
cance greater than 5. Once we identify the location of a 
source, we next determine whether the source is extended 
or unresolved, as seen by the SZA or the VLA, and 
whether this candidate is a single source, or a collection 
of nearby sources. Due to the complex sidelobe structure 
of the synthesized beam (see §2), nearby sources must be 
removed simultaneously from the interfcrometric data; 
we therefore fit any additional sources within 45" of the 

The Very Large Array is a facility of the National Radio As- 
tronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 



primary source location, roughly twice the synthesized 
beam width of the long baseline maps. 

Once we have determined all sources near the brightest 
in the map which are to be removed from the data, as well 
as their morphology (compact/extended), we solve for 
source properties by fitting to the multi-pointing visibil- 
ity data. For computational expediency, we describe the 
sources as functions with analytic Fourier transforms (see 
§4.2). The best-fit models are removed from the Fourier 
data, and the mosaics are regenerated. This process is 
repeated iteratively until there are no sources brighter 
than 5a in the significance maps. 

4.2. Models 

When fitting sources in the mosaicked maps, we model 
the intensity distribution of an unresolved source as a 
delta function, 

x) = lo d{x - Xo), (3) 

where u is the observation frequency, Iq is the intensity at 
frequency Uqi defined for the SZA as the center frequency 
of observations, and a is the spectral index. 

We model any extended source as a 2-dimensional 
Gaussian, i.e., up to a normalization, 

I{v)'X.Io{ — \ e ^'i' e ^'m^ , (4) 

where Ic and rUc are the coordinates of the centroid, 
is the orientation angle, and ai and am are the FWHM 
of the semi-major and semi-minor axes, respectively. In 
Eq. 4, the primed coordinate system {I', m') is related 
to the unprimed coordinate {l,m) via a rotation by the 
orientation angle 6. 

This model has an analytic Fourier transform provided 
we can take the primary beam of the SZA to be con- 
stant over the extent of the source, i.e., that the source 
is small compared to the primary beam. In practice, ex- 
tended models are fit only to pointings within 6' of the 
source, well within the FWHM of the primary beam. For 
a source 1' in extent, a factor of ~ 3 larger than the most 
extended source in our data set, this assumption leads to 
an error in the fitted fiux that is below 5%. 

4.3. Unresolved vs. Extended Sources 

Sources in the SZA maps are cross-checked against the 
VLA 5 GHz catalog, and a source is determined to have 
a counterpart at 5 GHz if a source exists within 8" of the 
SZA location, a small fraction of the highest resolution 
element of the SZA. If a counterpart is found, we consider 
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TABLE 2 
Survey Sensitivity 



Field 


Short Baselines 


Long Baselines 


All Baselines 


Area 


Name 


Minimum rms 


Median mis 


Minimum rms 


Median rms 


Minimum rms 


Median rms 


Covered 




(mjy/beam) 


(mJy/beam) 


(mJy/bcam) 


(m,]y/bcam) 


(mJy/bcam) 


(mJy/beam) 


(degree^) 


SZA4 


0.218 


0.264 


0.231 


0.283 


0.159 


0.193 


0.91 


DLS 


0.202 


0.238 


0.200 


0.251 


0.142 


0.173 


1.57 


NDWFS 


0.219 


0.239 


0.219 


0.249 


0.156 


0.172 


1.56 


SZA3 


0.213 


0.232 


0.218 


0.241 


0.153 


0.167 


1.74 


CMB fields 


0.152 


0.393 


0.153 


0.418 


0.108 


0.286 


1.91 



the source extended if its size at 5 GHz is greater than 
22.5"(the SZA long-basehne synthesized beam FWHM). 
In the rare event that the source has no 5 GHz counter- 
part (due presumably to source variability or resolution 
effects when observing extended sources) we check within 
8" of the SZA location in the 1.4 GHz NVSS or FIRST 
catalogs. When no counterpart is found in any of the 
catalogs (~ 3% of sources), we compare the flux at that 
location in the short and long baseline SZA maps to de- 
termine if it is extended. If the fluxes are consistent with 
each other the source is classifled as unresolved. 

4.4. Source Location Determination 

To minimize the number of degrees of freedom in each 

source fit, we fix the location of SZA sources with VLA 
counterparts to the VLA positions. We have verified that 
the locations of the VLA counterparts agree with loca- 
tions fitted to the SZA data to within 3", and indicate 
no systematic pointing offsets. 

Locations are not fixed for sources determined to be 
extended, and for very bright sources at 31 GHz. The 
dynamic range (ratio of the fitted source flux to the post- 
flt residuals) in the mosaicked maps is approximately 35, 
limited by the long-term pointing accuracy of the instru- 
ment and uncertainties in the primary beam. For sources 
brighter than 15 mJy, however, we can fit the locations 
accurately in the individual pointings, resulting in an 
overall dynamic range in the mosaics of 70, or a factor of 
2 improvement. We have verified that agreement in the 
fitted source location between pointings is at the arcsec- 
ond level. 

4.5. Spectral Index Fitting 

For source brighter than ^ 3.5 mJy. the SZA's large 
bandwidth allows for the simultaneous determination of 
the spectral index across 8 GHz. For fainter sources, 
however, we cannot meaningfully constrain the spectral 
index of the source from the SZA data alone. To reduce 
the number of degrees of freedom for these sources, we 
use the VLA 5 GHz observations (where available) as a 
second frequency, fixing the spectral index to 



-ln(/o//5GHz)/ln(;^o/5 GHz) 



(5) 



in Eq. 3. Note that these spectral indices are used only 
to reduce the residuals in the maps when extracting 
sources. Analysis of the spectral index distribution of 
these sources by comparison with low-frequency data is 
presented in §5.3. 

5. 31 GHZ POPULATION CHARACTERISTICS 
5.1. Source Sample 



Following the procedure outlined above, we identified 
209 sources at 31 GHz with snr greater than 5. From this 
sample, 40 of the sources were determined to be extended 
in the VLA data, yet none were seen to be resolved by the 
SZA, i.e. the FWHM of the fit Gaussians are all smaller 
than 22.5". In Table 1 of the Appendix we present the 
full list of sources. The brightest source we detect has a 
flux of 204 mJy. Counterparts were found in the VLA 
5 GHz catalog for 162 of the 171 sources with 5 GHz 
coverage, and 1.4 GHz NVSS counterparts were found 
for 157 of the 209 sources. All 1.4 GHz sources in the 
regions for which we have 5 GHz coverage were also found 
in the 5 GHz data. VLA 8 GHz coverage is available in 
the region for which we do not have 5 GHz coverage, 
and all 38 sources detected at 31 GHz in those regions 
have 8 GHz counterparts (Sharp et al. 2009). Of the 162 
sources with 5 GHz counterparts, 34 were found to have 
higher fluxes at 31 GHz than at 5 GHz. In addition, 8 
of the 157 sources with 1.4 GHz counterparts were found 
to be inverted. 

5.2. Source Count 

The differential source count per unit area as a function 
of flux, dN/dS, is typically described by a power law, i.e.. 



dN ^/S 



(6) 



where is a normalization parameter per unit area, 
5* is the source flux, and 7 is the power law index. In 
the analysis, we take Sq — 1 mJy and express the nor- 
malization as a number of sources per square degree. As 
models of radio source populations suggest a break in the 
count near 15 mJy, we fit a power-law only to sources be- 
low this limit, excluding 7 sources from the previous list. 
Subsequent analysis of the source count is performed on 
202 sources. 

5.2.1. Power law index, 7 

Given a set of Ns observed source fluxes {Si), we can 
solve for the underlying population parameters by max- 
imizing the likelihood of the dataset. 



Ns 

^p{Si\No,^,ai), 



(7) 



where Si is the observed flux of the source, and cTj is 
the map noise at the source location. In the presence of 
noise, an observed source flux S is the combination of a 
"true" flux St and noise N = S — St- The probability 
in Eq. 7 can therefore be written as an integral over all 
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possible pairs of St and N that will produce an observed 
flux Sf. 

p{Si\No,j,ai)<x / ps{St\No,^)pN{Si-St\ai)dSt 
Jo 

/•OO 

cx/ NoSr^e-('--^^^'/"^-dSt, (8) 

^0 

up to a constant of proportionality (e.g., Murdoch ct al. 
1973). It is clear that when p{Si\No, j, a^) is normalized, 
A^o drops out of the equation. Thus, this method can 
only be used to calculate the power law index 7; the rel- 
ative frequency of sources of different flux is related only 
to the shape parameter 7, and is independent of the total 
number of sources. In the noiseless case, it can be triv- 
ially shown that the resulting estimator for 7 meets the 
Fischer-Neyman criterion for a sufficient statistic, i.e., 
that it utilizes all the information about the shape pa- 
rameter contained in the data set (e.g., Crawford et al. 
1970). 

In the case of non-uniform noise, the above expression 
leads to a normalized likelihood that depends in a com- 
plicated way on the relative areas at different noise levels 
in the survey. However, because our sources are selected 
by snr and not flux, it is convenient to transform directly 
into the snr basis, leading to an expression for the like- 
lihood that is independent of the noise, and therefore of 
the relative areas at different noise levels: 

/•OO 

p(s,|7)(x / sr^e-(^--^*)'/2ds4, (9) 
Jo 

where s = S/a. We note in passing that while the kernel 
of Eq. 9 formally diverges as Si ^ 0, this is merely an 
artifact of the implicit assumption that the single-source 
hypothesis dominates the probability, i.e., that we are 
considering only sources well above the confusion limit. 
This assumption is consistent with the results quoted 
in §5.2.3, which indicate that the weakest source flux 
included in the analysis is approximately 14 times the 
confusion limit. In practice, we truncate Eq. 9 at an 
St which corresponds to a flux well above the confusion 
limit at the lowest map noise, and below which the kernel 
contributes negligibly to the integral. 

We estimate the uncertainty in the calculated value of 
7 by using a quadratic approximation of the likelihood 
about its maximum. We have verified through extensive 
simulation that this method produces an unbiased esti- 
mate of the power law parameter 7 and its associated 
uncertainty a^. 

5.2.2. Normalization, Nq 

We can invoke^ Bayes' theorem to write the likelihood 
of a set of parameters {A^o; 7} given the data D as 

p{No, j\D)(xp{D\No,-f) p{No, 7) 

cxp(7^|7Vo,7)p(iVo)p(7). (10) 

Assuming a uniform prior for Nq, p{Nq) = const, and 
integrating over 7, we have 

piNo\D) (X J piD\No,j)pij)dj. (11) 

The distribution ^(7) is just the likelihood derived in 
the last section, or for computational convenience, its 



Gaussian approximation. The probability of the data 
given Nq and 7 is the Poisson probability of the observed 
number of sources Ns, given the expected number of 
sources X{No,'y), 

p{D\No,,)c<^^, 

where X{No, 7) is derived by integrating dN/dS over the 
noise mask of the SZA survey (given in Table 2 of the Ap- 
pendix). The maximum likelihood estimate of the nor- 
malization is then just the value of A^o that maximizes 
p{No\D). 

5.2.3. Results 

Placing a significance cutoff on the source detections at 
5(7, the source count follows a power law with 7 = 2.18 ± 
0.12, and a normalization at 1 mJy of Nq = 32.1 ± 3.0 
deg~^ (in the 0.7— 15 mJy range). Integrating Equation 
6 with these parameters yields an integrated source count 
of 

N{>S) = (27.2 ± 2.5) deg-2 x iS„uv)-'-'^^°-'^ (12) 

As discussed in §2.5, the effective noise in the interior 
of the mosaicked maps is nearly uniform, but increases 
rapidly towards the edges. Because of this nonunifor- 
mity, the 5a significance threshold used to select sources 
independent of noise level does not translate simply into 
a uniform completeness limit across the full area of the 
survey. From the noise distribution, the sample is ex- 
pected to be 98% complete above 1.4 mJy over a 4.3 
square degree area, but 98% complete above 5 mJy over 
the full 7.7 square degree area. Here, we define the 98% 
completeness limit as the fiux above which the probabil- 
ity of missing a source due to noise fluctuations is < 2%. 

5.3. Spectral Index Distribution 

We calculate spectral indices for all sources in the sam- 
ple (including those brighter than 15 mJy) with counter- 
parts in the 5 GHz or 1.4 GHz VLA catalogs, where the 
spectral index a is defined as in Eq. 3. For each source, 
we construct the probability distribution for the spectral 
index by sampling from the flux distributions at each fre- 
quency. In the presence of noise, the conditional prob- 
ability of the "true" 31 GHz flux St, given an observed 
flux S is 

p{St\S,j,(T) oc psiSth)pN{S - St\a), (13) 

which is just the kernel of Eq. 8. For 7 > 0, Ps{St\j) 
increases as 5* — > 0, and the observed flux S is therefore 
generally "boosted" by the presence of noise. Marginal- 
izing this kernel over 7 using the likelihood derived in 
§5.2.1 yields the distribution of St- For the VLA fluxes, 
no correction for boosting has been made. This process 
is repeated for each source, and the resulting probabili- 
ties axe co-added to form the spectral index distribution 
for all sources with counterparts. 

By correlating with the NVSS catalog, as described 
in §4.3, we flnd 7 instances where the higher resolution 
of the SZA and VLA 5 GHz data, compared to NVSS, 
results in the detection of two sources associated with 
a single NVSS source. When calculating 1.4/31 GHz 
spectral indices, we treat these seven pairs as a single 



source. As a result, we calculate 1.4/31 GHz for 150 out 
of 202 sources, shown in the top panel of Figure 3. The 
distribution peaks at a spectral index of a 0.7, con- 
sistent with standard synchrotron emission, as expected 
for most sources bright at 1.4 GHz. We note that of 
the sources with 1.4 GHz counterparts, approximately 
5% are inverted, i.e., have a < 0. Accounting for sources 
undetected at 1.4 GHz will likely increase this fraction. 
For these sources, we use the NVSS completeness limit 
('^ 3.5 mjy) to construct upper limits for their spectral 
indices. Figure 3 illustrates the shift toward more in- 
verted spectral indices when these limits are included. 

We repeat the same calculation with the 162 sources 
with 5 GHz counterparts, noting that only 171 of the 
209 sources are covered by the 5 GHz follow-up. The re- 
sulting spectral index distribution for these 162 sources 
is shown in the bottom panel of Figure 3. The dis- 
tribution peaks at a value of a ~ 0.8. Of these 162 
sources, we find that 14% are inverted. Incorporating 
limits (S's GHz ^ 0-33 mJy) for nine sources which were 
not detected at 5 GHz has a minimal effect on the dis- 
tribution, shifting it marginally towards more inverted 
sources. 

We stress that information about the intrinsic spectral 
index distribution of this population of sources can only 
be inferred from the observed distributions by properly 
accounting for the selection thresholds at each frequency. 
For example, we can find particular choices of intrinsic 
spectral index distribution that reproduce both distribu- 
tions shown in Figure 3 (though in general we have no 
reason to expect a single distribution to describe both 
the 1.4/31 and 5/31 GHz spectral indices). We therefore 
refrain from drawing any conclusions about frequency 
evolution in spectral indices from these data. 

6. DISCUSSION 

In Figure 4, we plot the SZA dN/dS measurement over 
the range 0.7 — 15 mJy. In the top panel of the figure, 
we show results from 31 GHz-selected sources reported 
by OVRO/BIMA (Coble et al. 2007), CBI (Mason et al. 
2003), DASI (Kovac et al. 2002), and the VSA (Cleary 
et al. 2005). In the bottom panel we plot the projection 
from de Zotti et al. (2005), and from Mason et al. (2009), 
based on targeted 31-GHz foUowup of a 1.4 GHz-selected 
sample of sources from the NVSS catalog. 

6.1. Comparison with Cluster Source Measurements 

The only prior measurement of 31-GHz selected 
sources complete to a comparable flux level is the Coble 
et al. (2007) sample of '-^ 100 sources brighter than 1 mJy, 
found serendipitously in targeted observations of mas- 
sive clusters of galaxies. Since the vast majority of these 
sources lie within the inner regions of clusters, these 
measurements permit a direct comparison of the clus- 
ter population with the 31-GHz field-source population 
(this work). By integrating the field-source count from 
the SZA and the cluster source count from Coble et al. 
(2007) over the flux range 0.7— 15 mJy, we calculate the 
overdensity of radio sources within the central arcminute 
of clusters to be a factor 6.8 ± 4.4, consistent with their 
estimate. 

6.2. Comparison with Field Source Measurements 
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Spectral Index 

Fig. 3. — Top: Spectral index distribution for 31 GHz-selected 
sources detected with the SZA relative to their 1.4 GHz flux seen by 
NVSS. Red histograms denote sources with identified counterparts 
while blue histograms include upper limits for undetected sources, 
assuming the NVSS completeness limit of 3.5 mJy. Bottom: Same 
histograms but for counterparts found in the 5 GHz VLA follow-up 
data, with a limiting flux of 0.33 mJy at 5 GHz. 



From the top panel of Figure 4, we note the general 
agreement between the SZA result and prior 31-GHz 
measurements, at the higher completeness levels of those 
experiments (> 4 mJy). Agreement with the field-source 
count of Coble et al. (2007) is likewise good, although 
their sample in non-cluster fields consists of only four 
sources. From the bottom panel of Figure 4, we see that 
the SZA result is also consistent with the projections de 
Zotti et al. (2005) and Mason et al. (2009) (hereafter 
M09) for fiuxes greater than ~ 4 mJy. At low fiux, how- 
ever, they deviate significantly. 

The M09 projection is from a study of 3165 sources 
selected from the NVSS catalog at 1.4 GHz and reob- 
served at the same central observing frequency as the 
SZA. Adopting the 1.4 GHz source distribution of Hop- 
kins et al. (2003), and an intrinsic 1.4/31 GHz spec- 
tral index distribution inferred from their own obser- 
vations, they estimate the integrated 31-GHz source 
count over the same fiux range as the SZA results to 
be N{>S) = (16.7 ± 1.7) deg-2 x(S'„jy)-o-8"±o o7. This 
is inconsistent with the count we determine directly at 
31 GHz (§5.2.3) in both the normalization and the power 
law index, as shown in the bottom panel of Figure 4. 

As a consistency check, we generate simulated source 
populations under the set of assumptions outlined in 
M09, apply the source extraction algorithm detailed in 
§4.1 to the simulated data, and apply the formalism of 
§5.2 to calculate the expected N{>S) at 31 GHz. In 
particular, we generate a list of 1.4-GHz source fluxes 
over the range 25 /xJy — 1 Jy using the source count of 
Hopkins et al. (2003), and assigning spectral indices ac- 
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Fig. 4.— Top: Measurements of the 31 GHz dN/dS from this 
work and prior experimental results from OVRO/BIMA (Coble 
et al. 2007), CBI (Mason et al. 2003), the VSA (Cleary et al. 
2005), and DASI (Kovac et al. 2002). Bottom: Comparison of 
the SZA dN/dS to projections from lower frequencies by de Zotti 
et al. (2005) and Mason et al. (2009). 



cording to M09 to extrapolate these fluxes to 31 GHz. 
Sources are then assigned random locations on a noise 
map identical to the actual SZA survey coverage, with 
appropriate Gaussian noise added to the 31 GHz fluxes. 
We select only sources which would have been detected 
by the SZA, i.e., with snr greater than 5, and calculate 
the integrated source count. This procedure is repeated 
for 100 realizations of source populations, resulting in 
Ni>S) = (16.li^|)deg-2 x(5„^^)-0-86±o.i8^ This re- 
sult is consistent with the prediction of M09, and demon- 
strates that our methodology is robust to the experimen- 
tal details of the SZA survey, or differences in source 
selection between the two measurements. 

We have also investigated whether this discrepancy 
can be due to differences in resolution between the 
SZA and the VLA in its NVSS configuration. How- 
ever, of the 209 sources detected with the SZA, we 
note that in only 7 cases are multiple sources close 
enough to appear as a single source in the NVSS 
beam. Treating these cases as single sources, we ob- 
tain iV(>5') = (26.1±2.3)deg-^ x(S'„jj/)-i i9±o i3, an 
insignificant change from our nominal result (Eq. 12). 

To reconcile the two measurements, we postulate a 
change in the 1.4/31 GHz spectal index distribution for 
sources below the 1.4 GHz flux limit of the M09 study. 
Although M09 find no significant difference between the 
spectral index distributions of 1.4 GHz sources with flux 



> 10 mjy and sources with flux < 10 mJy, we note that 
their analysis is necessarily limited to sources brighter 
than the NVSS completeness limit of 3.5 mJy. Sources 
below this limit, however, constitute a large fraction of 
the 1.4 GHz source population used to predict the 31- 
GHz dN/dS, and small changes in the assumed spectral 
index distribution of this population can have a signif- 
icant impact on the source population at 31 GHz. We 
note also that the SZA data suggest this hypothesis ab 
initio; while the M09 distribution predicts that fewer 
than 10% of sources detected by the SZA at 31 GHz 
would lack 1.4 GHz counterparts > 3.5 mJy, fully 25% 
of SZA sources lack counterparts in the NVSS catalog. 

To test this hypothesis, we repeat the simulations 
described above, using the M09 spectral index distri- 
bution for sources above the NVSS completeness limit 
(3.5 mJy), and a separate spectral index distribu- 
tion for sources with 1.4 GHz fluxes below 3.5 mJy. 
For the dim sources, we assume a distribution whose 
shape is identical to that of M09, but shifted with re- 
spect to it. A good match to our data can be ob- 
tained by shifting the M09 spectral distribution by 
0.35 to more inverted spectra, resulting in a predicted 
N{>S) = (25.5 ± 2.9) deg-^ x (S'™j^)"i °6±o i4. This 
distribution has a peak spectral index of ^ 0.7, and pre- 
dicts that ^ 7% of dim sources are inverted between 1 .4 
and 31 GHz. This ad hoc model also accounts for the ob- 
served fraction of SZA sources that fall below the NVSS 
detection threshold. In addition, this model provides a 
better fit than the single-component model to the ob- 
served 1.4/31 GHz spectral index distribution shown in 
Figure 3. Note that while this ansatz agrees well with 
our data, a physically motivated model of faint source 
spectra is likely to be much more complex. 

FoUowup studies of faint 1.4-GHz sources at higher fre- 
quency have indeed found evidence for a fiattening of the 
spectral index distribution near the mJy level. Prandoni 
et al. (2006) find that for sources whose 1.4 GHz fiux is 
greater than 4 mJy, no sources are inverted between 1.4 
and 5 GHz, but that for sources dimmer than 4 mJy at 
1.4 GHz, roughly 10% of sources are inverted, compa- 
rable to the fraction of inverted sources we find between 
1.4 and 31 GHz. In addition, the study of Donnelly et al. 
(1987) conclude that the median 1.4/5 GHz spectral in- 
dex is roughly 0.75 down to a 1.4 GHz fiux of 0.25 mJy, 
and that the fraction of fiat spectrum sources increases 
from 22% at > 0.5 mJy to 41% in the 0.25 - 0.5 mJy 
range. 

7. CONCLUSIONS 

We present a sample of 31 GHz-selected sources 
from a 7.7 square degree survey obtained with the 
Sunyaev-Zel'dovich Array. We identify 209 sources 
at > 5(7 detection significance, ranging in fiux 
from 0.7 — 204 mJy. A maximum likelihood de- 
termination of the integrated source count results 
in N{> S) = (27.2 ± 2.5) deg'^ x {S^jy)-^-'^^^°-^^ in the 
fiux range 0.7 — 15 mJy. Comparison with a measure- 
ment of 31 GHz sources towards massive galaxy clusters 
leads to an overdensity of 6.8 ± 4.4, for sources within 
the central arcminute of massive clusters, relative to field 
sources. 

Of the existing source samples selected at 31 GHz, the 
SZA sample represents the only unbiased study of sources 
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valid to < 1 mjy. This sample will therefore be useful 
in refining source models that are currently constrained 
by measurements at much higher flux. Since sources 
are also the most significant foreground for small-scale 
temperature anisotropy experiments, these data will be 
useful for the accurate inference of cosmological infor- 
mation from the current generation of centimeter-wave 
CMB anisotropy and SZ cluster survey measurements. 
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APPENDIX 



TABLE 1 
SZA 30 GHz Sources 



(J2000) 



5 

(J2000) 



St 31 GHz 

(mJy) 



(mJy) 



Ss GHz Si. 4 GHz' 

(mJy) (mJy) 



09:16:37.50 29:50:35.53 14.60l^ J^ 
09:18:01.17 30:01:28.00 

09:21:30.36 30:06:01.78 6.43; 

09:20:28.42 29:47:57.00 4.38+ 



8 70+°-^^ 

°- '"-0.17 
+0.17 
0.18 



14.60 93.80 

8.70 55.13 

6.44 37.54 

4.39 30.12 



31.0 
32.0 
28.7 
19.0 



51.1 
93.7 
64.4 
10.2 



^ Flux corrected for boosting (i.e., St in Eq. 13). Errors do not reflect uncertainties in the absolute flux calibration 
^ Bcst-flt flux, not accounting for boosting (i.e., S in Eq. 13). 
^ Significance of detection in the composite mosaics 
^ VLA 5 GHz fiux. Map rrns is 68/^Jy/beam 

1.4 GHz Flux from NVSS catalog 
^ Single source in NVSS catalog. 31 GHz fluxes were combined in the 1.4/31 GHz spectral index analysis 
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TABLE 1 — Continued 







a a 
^t3l GHz 


b 

^31 GHz 




05 GHz 


G e 
01. 4 GHz 


(J2000) 


(.12000) 


(m,Jy) 


(mjy) 




(mJy) 


(mJy) 


09:22:35.52 


29:53:48.95 




4.09 


28.25 


4.1 


< 3.5 


09:18:33.50 


30:15:56.65 




3.66 


27.89 


17.1 


66.8* 


09:18:32.51 


30:16:09.89 


^■^■^-0-18 


2.48 


24.46 


9.2 


66.8* 


09:23:29.98 


30:11:07.68 


4-58-o'.20 


4.60 


23.68 


18.8 


35.5 


09:22:01.24 


30:14:11.86 


q yr + O.ie 
•5- 'O_o.l7 


3.77 


23.64 


3.1 


34.5 


09:21:05.11 


30:00:34.09 


3 44+0 15 


3.46 


22.59 


18.4 


52.8 


09:23:31.36 


29:49:27.80 


o 2Q+017 
•^■^"-0.18 


3.31 


17.85 


5.4 


48.4* 


09:23:31.88 


29:49:14.78 


1 r;4+0.17 


1.58 


19.84 


14.6 


48.4* 


09:20:37.26 


29:56:26.93 


9 cq + (^15 


2.55 


16.59 


5.8 


14.8 


09:18:11.73 


29:50:16.22 


q 1 9 + 0.15 
•^■1^-0.17 


3.14 


14.38 


15.3 


63.7 


09:16:17.99 


29:40:40.40 




6.14 


14.07 


12.3 


24.8 


09:19:00.02 


30:21:53.50 




2.45 


12.80 


0.9 


< 3.5 


09:17:58.65 


29:55:31.69 


'^-0.15 


1.77 


12.19 


4.6 


8.0 


09:19:11.46 


30:07:33.91 


1 qq+0.17 


1.97 


11.43 


4.2 


8.2 


09:15:32.97 


30:15:56.24 


9 -77+0.26 
^- ' ' -0-27 


2.83 


10.74 


4.2 


< 3.5 


09:16:40.53 


30:13:06.15 




1.88 


10.87 


0.8 


< 3.5 


09:19:52.47 


29:56:25.35 


-o.n 


1.12 


10.01 


4.6 


20.9 


09:19:40.43 


29:59:32.82 


1 48+016 


1.52 


9.99 


7.6 


25.2 


09:16:51.27 


29:41:24.43 


q q9 + 0.34 


3.39 


9.83 


3.1 


< 3.5 


09:18:11.97 


29:58:36.63 




1.40 


9.61 


1.5 


4.1 


09:16:25.35 


29:52:16.26 




1.15 


8.68 


3.2 


11.0 


09:21:27.65 


29:58:11.85 


1 46+°-'-'' 


1.49 


8.44 


3.8 


10.9 


09:23:39.94 


30:06:30.24 


1 74+0.22 
^- ' ^-0.23 


1.80 


8.54 


12.6 


39.7 


09:20:18.31 


29:48:10.68 


1 io+°-i^ 


1.15 


8.19 


0.4 


< 3.5 


09:20:41.00 


30:10:12.89 


T 97+0.16 
^■^ ' -0.18 


1.32 


8.15 


2.4 


4.1 


09:18:57.85 


30:21:07.24 


1 5(5+0.19 


1.61 


7.67 


5.6 


32.1 


09:18:38.44 


30:18:32.21 


1 9Q + 0.17 
-^■''''-0.18 


1.34 


8.01 


4.8 


11.6 


09:19:29.27 


29:59:57.34 


T 9q + 0.15 
-'-■'^•J-0.17 


1.28 


7.96 


8.2 


8.7 


09:21:04.44 


30:20:31.25 


T q9 + 0.17 
l-''^-0.18 


1.37 


8.00 


2.9 


4.2 


09:16:26.72 


29:59:23.77 


1 4r;+0-15 
l-^°-0.16 


1.49 


7.76 


1.2 


< 3.5 


09:17:01.04 


30:01:26.81 


1 f;7+0.16 
-'-■°'-0.17 


1.61 


7.63 


4.9 


16.5 


09:23:45.88 


29:43:09.36 


2 4S+0-34 


2.58 


6.77 


14.6 


47.9 


09:23:34.35 


30:12:05.13 


1 45+°-^° 

i-*O_0.21 


1.51 


7.09 


1.5 


< 3.5 


09:22:06.74 


30:28:58.82 


4.36l°:f, 


4.57 


6.79 


6.7 


3.5 


09:23:03.47 


29:59:37.21 


0.98t°:l^e 


1.03 


6.72 


< 0.33 


< 3.5 


09:20:50.45 


30:17:35.77 


1.03l°;l| 


1.09 


6.38 


4.6 


10.3 


09:16:49.95 


30:11:56.08 




1.00 


6.51 


< 0.33 


< 3.5 


09:15:49.44 


29:50:02.20 


0.94t°e:l? 


1.00 


6.11 


3.1 


5.9 


09:18:14.89 


30:13:44.34 


07+0.16 
u.y ( _o.i9 


1.03 


6.18 


4.9 


14.1 


09:23:19.44 


29:46:07.68 


97+o.is 

u.y 1 _Q 20 


1.05 


5.58 


0.6 


< 3.5 


09:22:16.60 


30:15:07.99 




1.10 


5.52 


1.0 


3.6 


09:22:01.25 


29:56:01.58 


'°-0.16 


0.81 


5.42 


1.0 


< 3.5 


09:17:06.56 


29:50:23.21 


n 7K+0.16 

'°-0.17 


0.84 


5.43 


4.5 


16.8 


21:34:12.54 


25:20:05.00 


100.74+°-|° 


100.74 


243.60 


145.0 


123.4 


21:28:19.99 


25:06:39.95 


13.6812:11 


13.68 


60.65 


55.1 


126.5 


21:27:19.43 


25:12:50.59 


4.421°;!; 


4.44 


39.38 


4.3 


68.8* 


21:27:18.69 


25:13:12.58 


, 04+016 


3.05 


34.28 


16.2 


68.8* 


21:27:18.44 


25:12:20.60 


1 67+°-" 

-0.17 


1.71 


24.01 


9.1 


2.1 


21:33:04.62 


25:20:47.99 


Oq + 0.18 
0-*3_(, jg 


6.85 


37.60 


15.7 


33.2 


21:28:52.38 


24:39:06.47 


8-16-0.24 


8.17 


35.33 


9.8 


4.8 


21:28:53.25 


25:20:21.47 


„ 41+0.17 


3.43 


19.31 


6.9 


11.1 


21:32:08.45 


25:12:09.99 


2.86l°:1^7 


2.88 


18.25 


3.6 


4.5 


21:26:39.91 


25:03:49.19 


2.37l°;l? 


2.39 


14.72 


4.2 


< 3.5 


21:33:57.20 


25:15:56.97 


9 09+0.21 
^-S^-0.21 


2.85 


13.96 


2.6 


4.4 


21:32:30.18 


25:02:11.37 


9 „.,+0.16 
' -0.17 


2.09 


12.57 


5.5 


15.1 


21:31:53.81 


24:39:33.05 


9 qn + 0.20 
z.o.^_Q 21 


2.43 


11.81 


6.7 


10.1 


21:30:44.02 


24:49:33.86 


9 (ic + O.lS 


2.08 


10.59 


7.8 


20.8 


21:29:57.88 


25:07:25.01 


1 cq+O.ie 
1-OJ_0.17 


1.57 


9.31 


4.0 


7.4 


21:27:21.49 


24:58:17.25 


■I rc: + 0.15 
i-O'>_0.17 


1.59 


9.47 


1.1 


< 3.5 


21:27:00.95 


24:52:26.77 


i.49i°„:i? 


1.53 


9.35 


1.6 


< 3.5 


21:33:55.85 


25:12:29.07 


1.92l°-2? 


1.97 


9.42 


3.3 


2.5 


21:33:36.18 


25:06:17.37 


1.38l°;l? 


1.43 


9.08 


8.6 


30.6 


21:27:55.18 


25:01:44.83 


, qo + 0.16 


1.42 


8.75 


1.4 


< 3.5 


21:33:10.62 


25:15:56.43 


1 34+0-16 

1-J1_0.17 


1.38 


8.12 


2.5 


5.0 


21:33:01.83 


24:53:25.04 


, or, + 0.16 
1--JU_q 17 


1.34 


8.09 


4.9 


3.1 


21:31:47.69 


25:18:27.96 


-1 q-i+0.16 


1.35 


8.05 


4.0 


9.1 
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a a 
^t3l GHz 


b 

^31 GHz 


snr" 


05 GHz 


G e 
01. 4 GHz 


(J2000) 


(J2000) 


(m,Jy) 


(mjy) 




(mJy) 


(mJy) 


21:27:12.18 


25:10:37.86 




1.22 


7.50 


10.9 


39.9 


21:32:52.46 


24:55:34.07 


1.53l°;l? 


1.56 


6.72 


2.7 


8.7 


21:28:09.99 


24:41:36.19 


1 OO + 0.18 


1.34 


7.06 


1.8 


< 3.5 


21:31:17.47 


24:42:33.05 


1 00+°-^'' 

l.UU_Q^g 


1.06 


6.78 


< 0.33 


< 3.5 


21:31:29.21 


24:38:37.85 




1.43 


6.64 


4.3 


11.0 


21:27:36.49 


24:38:47.79 


1 4Q+0-23 


1.56 


6.53 


1.9 


3.0 


21:28:51.49 


24:58:10.41 


-0.18 


1.03 


5.89 


1.7 


3.3 


21:28:44.81 


25:06:25.11 




1.05 


6.26 


1.7 


< 3.5 


21:28:13.46 


25:26:28.12 


-'-■'''^-0.36 


2.02 


6.29 


4.2 


11.7 


21:27:48.05 


25:14:38.07 


0.92l«;i« 


0.98 


5.73 


1.1 


< 3.5 


21:30:51.28 


24:51:39.28 


o.8o;°:i| 


0.87 


5.98 


11.8 


47.3 


21:33:26.51 


25:22:15.72 


u.,^z_Q 22 


1.02 


5.86 


7.5 


28.3 


21:28:13.38 


24:48:42.27 


gi+o-" 


0.98 


5.74 


3.8 


< 3.5 


21:27:08.59 


25:19:08.56 


0.96iS:l« 


1.03 


5.84 


< 0.33 


< 3.5 


21:29:54.40 


24:54:40.10 


0.83tS;l? 


0.89 


5.77 


5.8 


21.5 


21:27:22.57 


24:52:52.14 




1.01 


5.79 


< 0.33 


< 3.5 


21:30:30.76 


25:16:35.42 




0.91 


5.38 


4.4 


11.3 


21:29:06.54 


25:06:34.41 


U.»i_o.l7 


0.98 


5.41 


1.5 


4.2 


21:30:53.15 


25:15:48.59 


1 24+°i6 


1.28 


5.43 


1.3 


3.0 


21:33:15.39 


25:28:58.52 


q erq+0.74 


3.86 


5.20 


37.6 


39.1 


14:29:10.22 


35:29:47.05 




20.32 


104.59 


26.9 


24.4 


14:34:34.22 


35:10:09.53 


14.50l°;f3 


14.51 


47.02 


36.2 


69.6 


14:34:00.30 


35:11:52.30 




8.09 


45.83 


10.1 


12.2 


14:27:58.74 


34:59:19.92 


6 95+"- 

"■''•■'-0-18 


6.96 


39.29 


5.7 


6.8 


14:26:32.19 


35:08:14.88 




5.50 


31.69 


30.1 


95.0 


14:25:41.98 


34:58:39.37 




5.71 


23.23 


35.6 


154.7 


14:25:40.39 


34:58:05.78 


7 04+0-34 
' ■^^-0.34 


7.98 


20.96 


29.2 


67.4 


14:31:34.55 


35:15:11.30 


^■^^-0-17 


2.94 


20.24 


21.2 


76.0 


14:34:13.86 


34:55:05.49 


o 04+0. 18 

•^■''^-o.ig 


3.37 


18.35 


2.1 


4.8 


14:27:59.15 


34:55:21.07 


7Q+0.17 


3.80 


13.44 


7.6 


34.1 


14:33:15.26 


34:51:41.46 


q T 7+0.21 
•5-1 ' -0.23 


3.20 


16.26 


23.2 


90.8 


14:33:20.62 


34:50:34.85 


n oq + 0.25 
i-.00_Q 26 


1.90 


9.97 


18.9 


90.6 


14:27:02.24 


34:59:04.28 


^- '"-0.18 


2.78 


15.73 


8.2 


17.7 


14:27:41.14 


34:59:31.52 


9 cq+0.17 
'=-i'<5-0.17 


2.56 


14.61 


1.6 


< 3.5 


14:29:09.28 


35:03:25.16 


1 7n+0-i7 

^- '"-0.18 


1.74 


13.25 


7.8 


46.6* 


14:29:11.11 


35:03:20.78 


1 08+°-^'' 


1.14 


10.44 


2.7 


46.6* 


14:29:06.01 


35:11:15.69 


^■■^'J-o.ig 


2.33 


13.20 


6.2 


13.2 


14:28:50.46 


34:54:20.89 


^■^-'--O.lQ 


2.44 


12.98 


5.8 


9.5 


14:32:45.57 


34:55:38.10 


9 -1 -1 +0.17 
^■11-0.17 


2.14 


12.51 


0.6 


< 3.5 


14:33:09.89 


35:15:18.21 


1 09 + O-I7 
-'-■''^-0.17 


1.86 


11.63 


9.4 


29.2 


14:30:26.60 


35:19:20.71 


3 41+018 
■^■^-'--O-i? 


3.43 


11.03 


5.8 


20.8 


14:32:13.54 


35:09:41.03 


-1 YQ+O-ie 
1- '^-0.17 


1.76 


10.12 


6.6 


16.5 


14:29:22.90 


35:12:19.80 


1 70+0.17 

J^- ''5-0.18 


1.77 


10.48 


2.0 


< 3.5 


14:34:35.36 


35:26:22.24 


q n-i +0.31 
'^■'■'1-0-32 


3.08 


10.08 


14.3 


30.8 


14:27:17.93 


35:01:30.76 


1 75+"- 

•■'-0.18 


1.79 


9.98 


3.1 


4.9 


14:30:11.82 


35:00:20.16 


1 «n+0-i6 


1.83 


9.97 


5.4 


12.7 


14:26:51.48 


35:19:24.68 


■1 cq+0-16 

j^-oy-0.17 


1.62 


9.49 


5.7 


4.4 


14:29:53.81 


35:17:54.52 


1 92+0-113 
^•"^-0.17 


1.96 


9.36 


0.9 


< 3.5 


14:32:39.93 


35:11:58.74 


1 48+0-17 

l-^°-0.17 


1.53 


8.92 


0.8 


< 3.5 


14:32:37.98 


35:30:36.88 


1 70 + 0.20 
!■ ' ''-0-21 


1.84 


8.88 


7.0 


19.9 


14:28:50.71 


34:53:15.03 


-1 -..-) + (). 19 
1- '^-0.21 


1.76 


8.65 


2.1 


2.3 


14:34:21.94 


35:25:17.88 


1 64+0-20 

l'"^-0.22 


1.69 


8.48 


4.5 


9.1 


14:34:39.66 


35:08:27.84 


q 4q+0.38 
"•^•5-0.40 


3.52 


8.16 


2.8 


5.6 


14:34:34.17 


35:09:58.11 


9 99 + 0.30 
^■'^^-0.32 


2.31 


8.04 


36.2 


69.6 


14:31:12.34 


35:35:26.58 


4 19 + 0.50 


4.25 


7.57 


3.7 


6.8 


14:32:56.09 


35:33:39.50 


9 9C + O.32 
^■-^•■'-0.34 


2.35 


7.90 


17.5 


66.3 


14:27:26.98 


35:02:12.07 


1 26+0-1'' 
l-'^°-0.19 


1.31 


7.69 


0.6 


< 3.5 


14:28:51.10 


35:30:30.07 


1 39+0-20 


1.45 


7.24 


7.1 


23.1 


14:32:16.89 


35:02:48.92 


-1 -1-1+0-17 
^•1^-0.18 


1.17 


7.28 


0.7 


< 3.5 


14:28:25.47 


34:55:47.06 




1.13 


7.05 


4.8 


11.7 


14:29:23.66 


35:28:51.22 




1.32 


6.89 


2.9 


7.9 


14:26:34.41 


34:59:29.06 




0.94 


5.60 


3.6 


8.1 


14:34:32.42 


35:22:20.90 


. co + 0.26 
i.OO_Q 28 


1.63 


6.19 


< 0.33 


< 3.5 


14:30:09.94 


35:19:57.93 


o.9oi«:i| 


0.96 


6.02 


0.4 


< 3.5 


14:29:48.68 


35:17:47.91 


1 03+0-1'' 


1.09 


6.12 


5.2 


15.2 


14:29:05.66 


34:49:10.32 


l-o''-0.38 


1.98 


5.82 


11.7 


33.4 



12 



TABLE 1 — Continued 







o a 
^t31 GHz 


b 

^31 GHz 


snr'^ 


05 GHz 


e 
01.4 GHz 


(J2000) 


(.12000) 


(mJy) 


(m,Jy) 




(m,Iy) 


(mJy) 


14:33:01.39 


35:01:49.41 




0.89 


5.10 


0.7 


< 3.5 


14:32:59.91 


35:28:33.87 


0-85io.i9 


0.93 


5.41 


5.6 


17.8 


14:33:08.12 


35:31:50.30 


1 34+0.24 


1.43 


5.83 


1.3 


2.8 


14:31:00.42 


35:36:32.76 


. ot. + 0.74 
^■■'^-0.76 


5.57 


5.56 


12.8 


38.1 


14:32:39.56 


35:01:51.23 


0.85l«:l| 


0.92 


5.55 


6.0 


15.1 


14:28:34.05 


34:57:02.02 




0.91 


5.39 


< 0.33 


< 3.5 


14:34:03.67 


35:24:53.09 




0.89 


5.32 


< 0.33 


< 3.5 


14:29:34.77 


35:27:42.12 




0.87 


5.02 


0.6 


2.6 


02:18:42.20 


31:49:26.88 


1 47 47+0-22 
-^^ ' ' -0.22 


147.47 


662.28 


112.0 


56.8 


02:14:15.63 


32:03:50.44 




26.83 


131.45 


30.6 


11.7 


02:13:39.17 


31:59:27.36 


9 ^7+0. 16 


2.39 


23.85 


11.4 


62.8+ 


02:13:38.65 


31:59:07.72 


2 67+°-^'' 


2.69 


21.09 


2.2 


62.8+ 


02:13:38.33 


31:58:48.36 


^ O + 0.17 
^■-^^-0.17 


2.15 


23.64 


7.6 


< 3.5 


02:16:35.38 


32:00:22.55 


„ 90 + 0.15 


3.28 


21.08 


19.7 


47.5 


02:14:02.07 


32:06:13.47 


1 90+"- 

-^"'"-0-19 


1.94 


19.77 


0.9 


< 3.5 


02:14:02.15 


32:06:06.63 


1 Q1+O.I8 
-^■^-^-0-18 


1.95 


19.30 


0.6 


< 3.5 


02:16:45.58 


32:03:58.85 


3.2o;S:i« 


3.21 


16.59 


10.4 


19.5 


02:18:50.64 


32:03:28.96 


2.68t°:\l 


2.70 


16.43 


5.3 


9.8 


02:19:02.96 


31:59:47.82 


9 1 1 +0.16 
''•^^-0.17 


2.14 


10.50 


9.8 


27.5 


02:13:55.80 


31:52:59.74 




1.89 


10.36 


3.7 


10.0 


02:15:19.89 


31:52:42.50 


J-.DO_Q 20 


1.71 


9.53 


8.9 


19.1 


02:16:01.36 


31:48:22.72 


9 r:fi+0.29 

Z.OO_Q 3Q 


2.63 


9.23 


8.7 


12.7 


02:17:53.59 


31:50:14.93 


1 69+°-^° 


1.75 


8.84 


1.5 


< 3.5 


02:12:56.83 


32:02:22.54 




1.39 


8.41 


2.7 


5.7 


02:14:09.72 


32:06:39.97 


0.86t°o,ll 


0.95 


7.00 


3.9 


24.5 


02:13:31.80 


31:58:08.46 


i.i>J-0.17 


1.16 


6.84 


0.4 


< 3.5 


02:12:43.34 


31:58:56.93 




1.02 


5.76 


4.3 


9.6 


02:15:53.67 


32:03:03.79 




1.07 


6.21 


1.7 


3.0 


02:17:19.14 


31:58:11.28 


o.8ot"„:l« 


0.88 


5.42 


- 


< 3.5 


02:18:02.31 


32:02:59.07 


n 70+016 

^- ' •'-O-IS 


0.87 


5.11 


- 


< 3.5 


02:22:48.81 


33:57:19.00 


19 or + 0.13 
l^.>5O-0.15 


12.35 


87.84 


- 


189.3 


02:22:22.79 


34:00:15.65 


r co + O.Sl 
O.DO_Q 52 


5.74 


10.98 


- 


5.5 


02:22:46.21 


34:04:33.89 


, oq+0.31 


1.54 


5.18 


- 


15.1 


02:30:41.98 


33:56:21.46 




4.04 


20.28 


- 


35.3 


02:11:40.54 


33:04:17.31 


9 e:q+0.30 
z.o»_Q 32 


2.67 


8.56 


- 


21.4 


02:20:45.66 


32:57:05.64 


1 0+0. 24 
•'•J^ '5-0. 24 


3.18 


8.50 


- 


23.5 


02:24:36.11 


32:51:45.89 


4 43+0.30 


4.48 


14.09 


- 


8.5 


02:16:03.95 


33:03:45.08 


i.03t«:^23 


1.12 


5.27 


- 


< 3.5 


02:20:48.03 


32:41:06.53 


204.05i«:f3 


204.05 


417.06 


- 


921.6 


02:20:17.28 


32:38:53.09 




4.32 


16.28 


10.1 


10.2 


02:25:00.23 


32:31:36.95 


7 41+0.36 


7.45 


20.48 


- 


16.8 


02:24:40.90 


32:29:39.65 


1 91+0.27 
^■'^^-0.29 


1.33 


5.16 


- 


< 3.5 


02:24:27.44 


32:34:27.52 


o.83l°;l^ 


0.87 


4.72 


- 


2.3 


02:11:48.05 


32:09:32.41 




1.32 


5.05 


- 


< 3.5 


02:20:24.90 


32:10:45.30 


-^■^'-'-0.13 


1.46 


11.18 


7.3 


17.5 


02:20:11.45 


32:10:42.50 


1 4K+0.19 
^■*°-0.21 


1.52 


6.36 


5.9 


19.7 


02:24:40.70 


32:02:05.45 




6.24 


19.66 


- 


6.5 


02:24:40.56 


32:08:43.00 


1 9r + 0.13 
^■^°-0.13 


1.28 


9.94 


- 


7.1 


02:24:59.44 


32:14:04.83 


9 ,.+0.28 
^■'°*-0-29 


2.32 


6.30 


- 


25.3 


02:16:27.02 


32:08:15.51 


o.76t?,;l? 


0.81 


6.27 


0.7 


< 3.5 


02:16:28.14 


32:10:30.20 


n f'-) + "-12 

0.62_o.i4 


0.67 


5.48 


1.6 


3.4 


14:18:40.52 


35:30:02.51 


&.82t';,\t 


6.82 


47.72 


- 


91.5+ 


14:18:40.59 


35:30:23.21 




0.91 


5.88 


- 


91.5+ 


14:18:40.27 


35:29:41.86 


9 c:o + 0.21 
^..}0_Q 21 


2.57 


10.89 


- 


< 3.5 


14:18:40.86 


35:33:24.71 


T oQ + 0.18 


1.44 


8.22 


- 


27.4 


14:22:54.68 


35:32:08.97 


0.72iS;l* 


0.78 


5.28 


- 


3.4 


14:22:42.74 


35:31:48.44 


0-68to.l'7 


0.75 


5.05 


- 


4.5 


21:24:22.68 


24:55:12.95 


8.72iS:™ 


8.74 


29.65 


- 


4.5 


21:24:32.61 


25:03:17.39 


1 7t: + 0.32 
^' '^-0.34 


1.87 


5.31 




5.8 


14:22:38.77 


34:58:18.65 


, 07+0.24 
— 0.26 


1.47 


6.00 




< 3.5 


21:25:03.78 


25:25:12.85 




4.72 


8.56 




< 3.5 


21:24:34.23 


25:27:58.04 


U.»0_Q 20 


1.06 


4.41 




7.4 


02:11:43.10 


33:31:39.73 


^•'"-0.23 


4.29 


19.54 




12.6 


02:11:31.56 


33:23:02.49 


^■^^-0.30 


4.26 


14.80 




< 3.5 


02:19:30.26 


33:28:51.69 


57.66;°;f, 


57.67 


224.25 




100.7 


02:19:43.75 


33:30:28.66 


, ,7+0.20 
^•^ ' -0.22 


1.25 


6.10 




< 3.5 



TABLE 1 — Continued 



« |5 St3i GHz* S31 GHz'' snr" S5 gHz'' Si. 4 GHz" 

(.12000) (.12000) (m,Jy) (m,Jy) (mjy) (mjy) 



02:24:23.34 


33:32:30.01 


,, + 0.29 
•^"'•^-0.30 


3.29 


11.45 




18.3 


02:23:33.81 


33:24:35.41 


3 91+0 80 


4.11 


6.77 




54.5 


21:24:50.28 


26:01:50.36 


7 8>S+''-^3 

' ■°°~0.20 


7.86 


33.37 




8.6t 


21:24:50.34 


26:01:59.06 


1.00_Q 22 


1.91 


8.50 




8.6+ 


21:33:02.14 


26:04:17.99 


, ,„ + 0.29 
^••J'^_0.32 


1.46 


5.24 




< 3.5 


21:36:42.22 


25:57:41.47 


15.211°;™ 


15.24 


28.20 




10.4 


21:36:55.58 


25:56:43.63 


4.38l°:|? 


4.42 


14.43 




6.1 



TABLE 2 
SZA Area Coverage 



NoiM"> \'aluo y2 Area 
(m-l.\-) (dcg^) 



0.11 


0.000 


0.12 


0.014 


0.13 


0.037 


0.14 


0.080 


0.15 


0.219 


0.16 


0.638 


0.17 


2.578 


0.18 


3.629 


0.19 


4.012 


0.20 


4.281 


0.21 


4.504 


0.22 


4.697 


0.23 


4.868 


0.24 


5.023 


0.26 


5.323 


0.28 


5.554 


0.30 


5.758 


0.32 


5.940 


0.34 


6.106 


0.36 


6.257 


0.38 


6.396 


0.40 


6.525 


0.45 


6.807 


0.50 


7.036 


0.55 


7.226 


0.60 


7.385 


0.65 


7.514 


0.70 


7.613 


0.75 


7.684 



